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Abstract 

By implementing the new IR-improved Dokshitzer-Gribov-Lipatov-Ahareni-Parisi-Cahan- 
Symanzik (DGLAP-CS) kernels recently developed by one of us in the HERWIG6.5 envi- 
ronment we generate a new MC, HERWIRI1.0(31), for hadron-hadron scattering at high 
energies. We use MC data to illustrate the comparison between the parton shower gen- 
erated by the standard DGLAP-CS kernels and that generated by the new IR-improved 
DGLAP-CS kernels. The interface to MC@NLO, MC@NLO/HERWIRI, is illustrated. 
Comparisons with FNAL data and some discussion of possible implications for LHC phe- 
nomenology are also presented. 



1 Introduction 



In the era of the LHC, we must deal with requirements of precision QCD, which entails 
predictions for QCD processes at the total precision [1] tag of 1% or better, where by 
total precision of a theoretical prediction we mean the technical and physical precisions 
combined in quadrature or otherwise as appropriate. We accordingly need resummed 
C(a2^'"),0(a,aL"'),C(a2^"") corrections for n = 0,1,2, n' = 0,1,2, n" = 1,2, in the 
presence of parton showers, on an event-by-event basis, without double counting and with 
exact phase space. Essential large QED and EW effects [2H1] are handled by the simulta- 
neous resummation of large QED and QCD infrared(IR) effects, QED®QCD resummation 
[S] in the presence of parton showers, to be realized on an event-by-event basis by Monte 
Carlo (MC) event generator methods. Indeed, we know from Refs. [3],|1] that no precision 
prediction for a hard LHC process at the 1% level can be complete without taking the 
large EW corrections into account. 

In what follows, we present the first step in realizing our new MC event generator 
approach to precision LHC physics with amplitude-based QED®QCD resummation by 
introducing the attendant new parton shower MC for QCD that follows from our ap- 
proach. We recall that in Refs. [6] our resummed QED MC methods, based on the 
theory in Ref. [7], are already well developed and checked in LEPl and LEP2 precision 
physics applications. This means that what we do here will set the stage for the complete 
implementation, via MC methods, the QED(8)QCD resummed theory in which all IR sin- 
gularities are canceled to all orders in and a. As we will show directly, this new parton 
shower MC, which is developed in the HERWIG6.5 [8] environment and which we have 
called HERWIRI1.0(31) [9)^1 . already shows improvement in comparison with the FNAL 
soft px data on single Z production as we quantify below. On the theoretical side, while 
the explicit IR cut-offs in the HERWIG6.5 environment will not be removed here, our new 
shower MC only involves integrable distributions for its real emission so that in principle 
these cut-offs could be removed. We discuss this point further below as well. 

Our discussion here proceeds as follows. We first review our approach to resummation 
and its relationship to those in Refs. [inilll]. Section 3 contains a presentation of the 
attendant new IR-improved DGLAP-CS [T2l[T3] theory p3l[T5] . Section 4 features the 

^We stress that we arc completely replacing all the parton shower kernels in HERWIG6.5 that generate 
real QCD radiation and all the Sudakov form factors in HERWIG6.5 that realize the attendant virtual 
corrections with the new forms that follow from the results in Sections 2 and 3, together with the auxiliary 
functions required for these new forms, so that the parton shower physics, distributions and MC behavior 
are all fundamentally different from what is in HERWIG6.5. The attendant implementation has been 
carried out in agreement with and in close collaboration with B. Webber and M. Seymour, principal 
authors of the HERWIG6.5, who have also instructed us on the proper naming and references for the 
resulting program as we further highlight in the discussion below. 
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implementation of the new IR-improved kernels in the framework of HERWIG6.5 [S] to 
arrive at the new, IR-improved parton shower MC HERWIRI1.0(31). We illustrate the 
effects of the IR-improvement first with the generic 2— )'2 processes at LHC energies and 
then with the specific single Z production process at LHC energies. We compare with 
recent data from FNAL to make direct contact with observation. Section [S] summarizes 
our discussion. 



For reference purposes and to put the discussion in the proper perspective with regard 
to what has already been achieved in the relevant literature, we note that the authors 
in Refs. [T6l[T7] have argued that the current state-of-the-art theoretical precision tag on 
single Z production at the LHC is (4.1 ± 0.3)% = (1.51 ± 0.75)%(QCD) © 3.79{PDF) © 
0.38 ± 0.26{EW)% and that the analogous estimate for single W production is ~ 5.7%. 
We continue to emphasize that these estimates show how much more work is still needed 
to achieve the desired 1.0% total precision tag on these two processes, for example. 



2 QED(g)QCD Resummation 

We follow here the discussion in Refs. [5l[Tll[T5], wherein we have derived the following 
expression for the hard cross sections in the SM SU2L x Ui x SU^ EW-QCD theory 

suMir(qced) 1 f d^P2 d^<l2 -TT d^kj^ -p-r d?k' 



^Z^^^l^iy po qo 11 j.^^ 11 



X 



y Jy-{pi+qi-P2-q2-y.kij -y,k' i^)+DocED q (U^ h ■ h' h' 



where the new YFS-style [7] residuals /3n^m{ki, . . . , kn, k[, . . . , k'^) have n hard gluons 
and m hard photons and we show the final state with two hard final partons with mo- 
menta P2, q2 specified for a generic 2f final state for definiteness. The infrared functions 
SUMir(QCED), -Dqced are defined in Refs. [51114( 115]. Eq. ([1]) is an exact implementation 
of amplitude-based simultaneous resummation of QED and QCD large IR effects valid 
to all orders in a and in Og- When restricted to its QED aspect, it is the basis of the 
well established YFS MC approach [B] to precision multiple photon radiative effects that 
is well tested already in LEPl and LEP2 precision physics applications. Thus what we 
present in this paper, the first realization of the new parton shower MC for QCD that 
follows from the QCD aspect of ([1]), opens the way to the full MC implementation of 
all aspects of our QED©QCD resummatiom theory approach to precision LHC physics 
predictions. 
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The approach to QCD resummation contained in ([T]) is fully consistent with that 
of Refs. [inilTT] as follows. First, Ref. [18] has shown that the latter two approaches 
are equivalent. We show in Refs. [HI [15] that our approach is consistent with that of 
Refs. [To] by exhibiting the transformation prescription from the resummation formula 
for the theory in Refs. [TD| for the generic 2 — t- n parton process as given in Ref. [T2j 
to our theory as given for QCD by restricting Eq. ([T]) to its QCD component, where a 
key point is to use the color-spin density matrix formulation of our residuals to capture 
the respective full quantum mechanical color-spin correlations in the results in Ref. [T9] . 
For completeness, we let us recapitulate the essence of the attendant discussion here, as 
the arguments are not generally well-known. More precisely, to illustrate the relationship 
between our approach and that in Refs. [TU], we use as a vehicle Ref. [H], which treats 
the 2 — )■ n parton process in the resummation theory of Refs [Tn], working in the IR 
and coUinear regime to exact 2-loop order. The authors in Ref. [19] have arrived at the 
following representation for the amplitude for a general 2 — n parton process [f] at hard 

scale Q, /i(pi,ri) + f2iP2,r2) fsiPs^r^^) + fi{pi,ri) H h /n+2(Pn+2, ^^^+2), where the 

PiiTi label 4-momenta and color indices respectively, with all parton masses set to zero 
( so in our approach, we should have in mind that the masses of the quarks (see the 
discussion below) and the IR regulator mass of the gluon would all be taken to zero or, we 
could, as it is done Ref. [19], just set all masses to zero at the outset and use dimensional 
regularization to define both collinear and IR singular integrals) 

' c (2) 
= JWX:5l///F(cl){M, 

L 

where repeated indices are summed, and the functions J^^\ Sli, and H^^ are respectively 
the jet function, the soft function which describes the exchange of soft gluons between 
the external lines, and the hard coefficient function. The latter functions' infrared and 
collinear poles have been calculated to 2-loop order in Refs. [19]. How do these results 
relate to eq.(IT])? 

To make contact between eqs. (fT]l2|) . identify in the specific application Q'Q — )• Q"'Q" + 
m(G) in (HD A = Q, ^ = Q', h = Q", U = Q'" . {/s, ■ ■ ■ , fn+2] = {d, ■ ■ ■ , G™},in ©, 
where we use the obvious notation for the gluons here. This means that n = m + 2. Then, 
to use eq.© in eq.([T]), one simply has to observe the following: 

I. By its definition in eq.(2.23) of Ref. [IH], the anomalous dimension of the matrix Sn 
does not contain any of the diagonal effects described by our infrared functions 
SUMir(QCD) and Dqcd, where 

SUMir(QCD) = 2as'^BQCD + 2a,EQCD(i^max), 
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/d k ~ 

Dqcd = J ^S^CDik) [e-''-' - ^(i^^ax - k)] , (3) 

, where the real IR emission function 5'qcd(^) and the virtual IR function ^Bqcd 
are defined eqs. (77,73) in Ref. [H]. Note that ([1]) is independent of Kmax- 

II. By its definition in eqs. (2.5) and (2.7) of Ref. [19], the jet function jt^l contains 
the exponential of the virtual infrared function ag^BqcD, so that we have to take 
care that we do not double count when we use ([2]) in ([1]) and in the equations in 
Refs. [5l[Tll[T5] that lead thereto. 

When we observe these two latter points, we get the following realization of our approach 
using the results in Ref. [19]: In our result in eq.(75) in Ref. [H] for the contribution to 
([T]) of m-hard gluons for the process under study here, 

" — — / 11 + X2)i/2 ^^P^ + qi-P2-q2 -2^h) 

j=\ ' 1=1 

'{pi,qi,P2,q2,ki,- ■ ■ ,km) o~0~' ^'^) 



pm 



we can identify the residual p("^) as follows: 
P^'^\pi,qi,P2,q2,h,--- ,k^) = Y. , l-^li'-'' 

^— ^colors. smn i'> 



color s, spin 
spins, {ri},{r'^} L=IL'=1 



(5) 



where here we defined .1^^^ = e~"''^'^'3CD g^^id we introduced the color-spin density 
matrix for the initial state, i)'^'^, so that ^fr,}{r'} ~ ^{^ r2}{r' r'}' suppressing the spin 
indices, i.e., i)'^^ only depends on the initial state colors and has the obvious normalization 
implied by (jl]). Proceeding then according to the steps in Ref. [H] leading from (jl]) to ([1]) 
restricted to QCD, we get the corresponding implementation of the results in Ref. [T^ in 
our approach, without any double counting of effects. 



3 IR-Improved DGLAP-CS Theory 

We show in Refs. [TUdn] that the result Eq. ([T]) restricted to QCD allows us to improve 
in the IR regime the kernels in DGLAP-CS [T^ [T3] theory as follows, using a standard 
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notation: 



is 1 + 



z 



Gq 



Z A. Z 

P;^\z) = Fy,s{lG)e^'-\{z\l - zyo + (1 _ ,)2,7c|^ (6) 

where the superscript "exp" indicates that the kernel has been resummed as predicted by 
Eq. ([T]) when it is restricted to QCD alone and where 

as. 4:Cf . 7g , asCpfTT^ l^ 



fqilqj 



71 (3o ^2 TT ' 3 2' 

2 2 1 



7g 7g + i Iq 



1 7. + 2' 



a, 4Cg 7g , OsCg^vt^ 1 

7g = Cg— t = -^, = ^ + — — - 7^)' 

71 Pq Z 71 O Z 

= e-^'5<= + + (7) 

^^^^^^ 6CgFyfs{1g) 7g(1 + 7g)(2 + 7g) (1 + 7g)(2 + 7g) 

1 1 

+ ttt:^ TT^ r + 



2(3 + 7g) (4 + 7g) (2 + 7g) (3 + 7g) (4 + 7g) 
Fyfs{i)= ^ V ^7 = 0.57721566..., (8) 

r(i + 7) 

where T{w) is Euler's gamma function and C is Euler's constant. We use a one-loop 
formula for as{Q), so that 

2 

/3o = 11 - 

where n/ is the number of active quark flavors and Cp = 4/3 and Cq = 3 are the respective 
quadratic Casimir invariants for the quark and gluon color representations. 

For the sake of completeness, let us illustrate how one applies the result in ([T]) to 
obtain the results in ([6]). We use the example of Pgg for definiteness. We apply the QCD 
exponentiation master formula embedded in eq.([T]) to the gluon emission transition that 
corresponds to Pqq{z), i.e., to the squared amplitude for q — )■ q{z) + G{l — z) so that in the 
specialized case already discussed above one replaces everywhere the squared amplitudes 
for the Q'Q — > Q"'Q" processes with those for the former one plus its nG analogs with the 
attendant changes in the phase space and kinematics dictated by the standard methods; 
this implies that in eq.(53) of the first paper in Ref. [12] we have from the application of 
the QCD aspect of eq.([T]) the replacement (see Fig. [1]) 
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Pi P2il-z) Pa P2il-z) 

(a) (b) 
Figure 1: In (a), we show the usual process q — )■ g(l — z) + G{z)] in (b), we show its 
multiple gluon improvement q — )■ g(l — z) + Gi{^i) + ■ • • + Gn{C,n), z = J2j 0- 



1 n \^ 



^BA = PL^^Z{1-^)Y1 ^ 

spins -L 



B+C\ 



|2 



(9) 



|2 



spins 

where A = q, B = G, C = q and Va^b+c is the lowest order amplitude for q — 
G{z) + g(l — z), so that we get the un-normalized exponentiated result 

P,,{z) = C^Fy^<,(7^)el^'i±^(l - zV". (10) 

1 — 2 

We see immediately that the exponentiation has removed the unintegrable IR divergence 
at z = 1. For reference, we note that we have in ( ITOl) resummed the term^ (9(ln^(l — 
z)t^ag), n > i > k, which originate in the IR regime and which exponentiate. The 
important point is that we have not dropped outright the terms that do not exponentiate 
but have organized them into the residuals f3m in the analog of eq.([I]). 



The application of eq.([T]) to obtain eq. (fTO|) proceeds as follows. First, the exponent 
in the exponential factor in front of the expression on the RHS of eq.(IT]) when restricted 

■^Following the standard LEP Yellow Book [20] convention, we do not include the order of the first 
nonzero term in counting the order of its higher order corrections. 
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to QCD is readily seen to be , using the known results for the respective real and virtual 
infrared functions from Refs. [HlfTS], 



SUMjniQCD) = 2a,5R5QCD + 2asBQCD{K^ 



^ (nry "^^1 ^^^^ , ^ "''^ , (^sCf ,T^'^ l^^ (H) 

— ( ZLjp — tin — — hOi? — 1 + 



2 V E 27r TT ' 3 2- 

where on the RHS of the last result we have already applied the DGLAP-CS synthesization 
procedure as prescribed in Refs. p!^j[T5] to remove the collinear singularities, In Aq(^^/"^^ 



m 



q 



1, in accordance with the standard QCD factorization theorems |2T]. This means that, 
identifying the LHS of eq.(II]) as the sum over final states and average over initial states 
of the respective process divided by the incident flux and replacing that incident flux 
by the respective initial state density according to the standard methods for the process 
q — > q{l — z) + G{z), occurring in the context of a hard scattering at scale Q as it 
is for eq.(53) in the first paper in Ref. [12], the soft gluon effects for energy fraction 
< z = Kmax/E give the result, from eq.(IT]) restricted to QCD, that, working through to 
the /^-level and using q2 to represent the momentum conservation via the other degrees 
of freedom for the attendant hard process, 

^'^^^ P^Adtdz - cSUMiR.(QCD)(2) /"rj f A \iy.{p^-p^)+ f^^Oi^-- !^Socu{k)[e-'y^~l]} 

27r J J (27r)4 



+ 



(fp2(Pq2 

' ' J 

^ SUMir(QCD)(z) /" r g /"°° {iy{E^-E2)+t^'^^-- !^SciCJ^{k)\e-'y>'-l\} 

y^^V_(2vr)' 

+ ffhp^fki) r cfa-(-gi-^2-fc?)+r<^'"-4^5oCD(fc)[e"'^-''-l]} 

J h J-oo (27r) 

(12) 

where we set Ei = p^, 2 = 1,2 and the real infrared function S'qcd(^) is known as well: 



( IL 


\ 2 
- — I 1 


\kpi 


kp2) 



SQCDik) — — ^ I j I DGLAP-CS synthesized (13) 



and we indicate as above that the DGLAP-CS synthesization procedure as prescribed in 
Refs. [111115] is to be applied to its evaluation to remove its collinear singularities; we are 
using the kinematics of the first paper in Ref. [12] in their computation of Pba{z) in their 
eq.(53), so that the relevant value of is indeed Q^. It means that the computation 
can also be seen to correspond to computing the IR function for the standard t-channel 
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kinematics and taking | of the result to match the single line emission in Pcq- The two 
integrals needed in (fT2|) were already studied in Ref. [7]: 



iYFs{zE,Q)= / 

YFSVlq 



iYFsizE,ki) = / — e 

oo 



^ [i3/(^E-fcl)+/'=<^^ 4^5QCD(fc)(e-'^'=-l)] 



oo , . „ _ (14) 



'zE - ki 



)'-^''Iyfs{zE,0) 



When we introduce the results in f lT4|) into f|T2|) we can identify the factor 

/ (^"^ + / **"'"'^'(^''(;^)*-^'^) ll - / "'"-i^^'^^'^ + ^("2)- 

(15) 

where is the unexponentiated result in the first line of This leads us finally to 
the exponentiated result in the second line of by elementary differentiation: 

Pba = P°BAz'''FYFs{%)e-^''' (16) 



Let us stress the following. In this paper, we have retained for pedagogical reasons 
the dominant terms in the resummation which we use for the kernels. The result in the 
first line of ( IT2l) is exact and can be used to include all higher order resummation effects 
systematically as desired. Moreover, we have taken a one-loop representation of for 
illustration and have set it to a fixed- value on the RHS of f|T2|) . so that, thereby, we are 
dropping further possible sub-leading higher order effects, again for reasons of pedagogy. 
It is straight forward to include these effects as well - see Refs. [T^ITS] for more discussion 
on this point. Repeating the exhibited resummation calculation for the other kernels leads 
to the results in (|6]). The latter results have now been implemented by MC methods, as 
we exhibit in the following sections. 



We stress that the improvement in Eq. should be distinguished from the also 
important resummation in parton density evolution for the "2; — )■ 0" regime, where Regge 
asymptotics obtain - see for example Ref. [2^123]. This latter improvement must also be 
taken into account for precision LHC predictions. 



Let us now recall that already a number of illustrative results and implications of the 
new kernels have been presented in Refs. [ri l [T5 | l2^ which we summarize here as follows 
for the sake of completeness. Firstly, we note that the connection to the higher order 
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kernels in Refs. has been made in Ref. [2]. This opens the way for the systematic 
improvement of the results presented herein. Secondly, in the NS case, we find [13] that 
the n = 2 moment is modified by ~ 5% when evolved with Eq. from 2GeV to lOOGeV 
with Uf = 5 and Aqcd — 0.2GeV, for illustration. This effect is thus relevant to the 
expected precision of the HERA final data analysis [JSj. Thirdly, we have been able to 
use Eq. ([T]) to resolve the violation [271I2B] of Bloch-Nordsieck cancellation in ISR(initial 
state radiation) at 0{al) for massive quarks [29]. This opens the way to include realistic 
quark masses as we introduce the higher order EW corrections in the presence of higher 
order QCD corrections - note that the radiation probability in QED at the hard scale 
Q involves the logarithm ln{Q'^/m'^), and it will not do to set rriq = to analyze these 
effects in a fully exclusive, differential event-by-event calculation of the type that we 
are constructing. Fourthly, the threshold resummation implied by Eq. ([T]) for single Z 
production at LHC shows a 0.3% QED effect and agrees with known exact results in QCD 
- see Refs. [5l[30l[3T]. Fifthly, we have a new scheme [15] for precision LHC theory: in an 
obvious notation, 

where the primed quantities are associated with Eq. (E]) in the standard QCD factorization 
calculus. Sixthly, we have [S] an attendant shower/ME matching scheme, wherein, for 
example, in combining Eq. ([1]) with HERWIG [8], PYTHIA [32], MC@NLO [33] or new 
shower MC's [33], we may use either pr-matching or shower-subtracted residuals 

{f^n,m{ki, . . . , kn] k[, . . . , k'^)} to Create a paradigm without double counting that can be 
systematically improved order- by order in perturbation theory - see Refs. [5]. 

The stage is set for the full MC implementation of our QED(8)QCD resummation 
approach. We turn next to an important initial stage of this implementation - that of 
the kernels in Eq. (jSj). 

4 MC Realization of IR-Improved DGLAP-CS The- 
ory 

In this section we describe the implementation of the new IR-improved kernels in the 
HERWIG6.5 environment, which results in a new MC, which we denote by HERWIRIl.O, 
which stands for "high energy radiation with IR improvement" [35] 

Specifically, our approach can be summarized as follows. We modify the kernels in 
the HERWIG6.5 module HWBRAN and in the attendant related modules [36] with the 
following substitutions: 

DGLAP-CS Pab IR-I DGLAP-CS P'/i (18) 
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while leaving the hard processes alone for the moment. We have in progress [22] the 
inclusion of YFS synthesized electroweak modules from Refs. [6]for HERWIG6.5, HER- 
WIG++ [38] hard processes, as the CTEQ [39] and MRST(MSTW) [IQ] best (after 2007) 
parton densities do not include the precision electroweak higher order corrections that 
do enter in a 1% precison tag budget for processes such as single heavy gauge boson 
production in the LHC environment [3]. 



For definiteness, let us illustrate the implementation by an example [111112], which 
for pedagogical reasons we will take as a simple leading log shower component with a 
virtuality evolution variable, with the understanding that in HERWIG6.5 the shower 
development is angle ordered [H] so that the evolution variable is actually ~ E9 where 9 
is the opening angle of the shower as defined in Ref. [H] for a parton initial energy E. In 
this pedagogical example, which we take from Ref. [H], the probability that no branching 
occurs above virtuality cutoff Qo is Aa(Q^, Qq) so that 



which implies 



Aa(g',Q^)=exp 



«^ dt 



dz^Pbaiz) 
/TT 



(19) 



(20) 



The attendant non-branching probability appearing in the evolution equation is 



A(g^t) 



t = the virtuality of gluon a. 



The respective virtuality of parton a is then generated with 

A,{Q^t)=R, 

where -R is a random number uniformly distributed in [0, 1] . With (note /3o 
here, where Uc is the number of colors) 



2tt 



(21) 



(22) 



(23) 



we get for example 



dz— Pgciz) 



An 



2nhQ In 
2 1 
In 



1 

— ^ / dz- \ 
91] Jo 2 L 



z' + (1 - zY] 



(24) 
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so that the subsequent integration over dt yields 



1 dt 



Q2 3 t 60 In (jj) 



36o 
2 



In In 



In 



1 1( 


r 






f\n 




I In 


(Ql 

\ A2 



(25) 



Finally, introducing / into Eq. f l20|) yields 



exp 



36n 



ln(^ 



ln(^ 



ln(|f 



3bo 



If we now let Aa(Q2, t) = i?, then 



In 



A2 



2 

36o 



R 



which implies 



t = K 



2iQ 



2\ R'^ 



A2 



Recall in HERWIG6.5 [8] we have 



(26) 



(27) 



(2^ 
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-(llrie-lO), nf = 5 
2 

-BETAF 

3 



(29) 



where in the last line we used the notation in HERWIG6.5. The momentum available 
after a qq split in HERWIG6.5 [8] is given by 



(QLST \ 
_I 
QCDL3 J 



JIJBETAF 



(30) 
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in complete agreement with Eq. fl25]) when we note the identifications t = QQBAR^, A = 
QCDL3, Q = QLST. 

The leading log exercise leads to the same algebraic relationship that HERWIG6.5 
has between QQBAR and QLST but we stress that in HERWIG6.5 these quantities are the 
angle-ordered counterparts of the virtualities we used in our example, so that the shower 
is angle-ordered. 



Let us now repeat the above calculation for the IR-Improved kernels in Eq. ([6]). We 



have 



so that 



z'{i-zy^ + {i-zyz 



4:FYFs{lG)e' 







2^ ■ &oln(^j(7G + l)(7G + 2)(7G + 3) 

This leads to the following integral over dt: 

dt 4FY^5(7G)e^^/^ 
Q2 t 6oln(^)(7G + l)(7G + 2)(7G + 3) 



(31) 
(32) 



4:FYFs{lG)e 



7g/4 



Ei 1 



8.369604402 \ 

bo (7G + 1) (7G + 2) (7G + 3) bo\n{^) ) 

We finally get the IR-improved formula 

Aa{Q\t)=exp [- {F {Q')-Fit))], 

where 

4FYi.5(7G)e^°/^ 



F{Q' 



bo (7G + 1) (7G + 2) (7G + 3) 



Ei 1 



8.369604402 \ 



bn In 



mi 



(33) 



(34) 



(35) 



and Ei is the exponential integral function. In Fig. [2] we show the difference between the 
two results for Aa{Q'^,t). We see that they agree within a few percent except for the 
softer values of t, as expected. We look forward to determining definitively whether the 
experimental data prefer one over the other. This detailed study will appear elsewhere [13] 
but we begin the discussion below with a view on recent FNAL data. Again, we note 
that the comparison in Fig. [2] is carried out at the leading log virtuality level, but the 
sub-leading effects suppressed in this discussion will not change our general conclusions 
drawn therefrom. 

For further illustration, we note that for the q — t- qG branching process in HER- 
WIG6.5 we have therein the implementation of the usual DGLAP-CS kernel as fol- 
lows: 



12 



WMIN = MIN(ZMIN*(1. -ZMIN) , ZMAX* (1 . -ZMAX) ) 
ETEST = (1. + ZMAX**2) * HWUALF (5-SUD0RD*2 , QNOW*WMIN) 
ZRAT = ZMAX/ZMIM 
30 Zl = ZMIN * ZRAT**HWRGEN(0) 
Z2 = 1. - Zl 

PGQW = (1. + Z2*Z2) ( 
ZTEST = PGQW * HWUALF (5-SUD0RD*2, QN0W*Z1*Z2) 
IF (ZTEST .LT. ETEST*HWRGEN(1) ) GOTO 30 

where the branching of g to G at 2; =Z1 occurs in the interval from ZMIN to ZMAX set by 
the inputs to the program and the current value of the virtuality QNOW, HWUALF is the re- 
spective function for in the program and HWRGEN(J) are uniformly distributed random 
numbers on the interval from to 1. It is seen that Eq. fl5^ is a standard MC realization 
of the unexponentiated DGLAP-CS kernel via 

as{Qz{l - z))Pg,{z) = a,(Qz(l - z))l±il:iil! (37) 

where the normalization is set by the usual conservation of probability. To realize this with 

the IR-improved kernel, we make the replacement of the code in Eq. ( l36i) with the lines 
NUMFLAV = 5 

BO = 11. - 2 . /3 . *NUMFLAV 
L = 16./(3.*B0) 

DELTAQ = L/2 + HWUALF (5-SUD0RD*2 , QNOW*WMIN) *1 . 184056810 
ETEST = (1. + ZMAX**2) * HWUALF (5-SUD0RD*2 , QNOW*WMIN) 

* EXP(0.5*DELTAQ) * FYFSQ (NUMFLAV- 1) * ZMAX**L 
ZRAT = ZMAX/ZMIN 
30 Zl = ZMIN * ZRAT**HWRGEN(0) 
Z2 = 1. - Zl 

DELTAQ = L/2 + HWUALF (5-SUD0RD*2 , QN0W*Z1*Z2) *1 . 184056810 
PGQW = (1. + Z2*Z2) * EXP(0.5*DELTAQ) * FYFSQ (NUMFLAV- 1) 
* Z1**L 

ZTEST = PGQW * HWUALF (5-SUD0RD*2, QN0W*Z1*Z2) 
IF (ZTEST .LT. ETEST*HWRGEN(1) ) GOTO 30 

so that with the identifications 7^ = L, 5q = DELTAQ, FYFsilq) = FYFSQ(NUMFLAV - 1), 
we see that Eq. f p8|l realizes the IR-improved DGLAP-CS kernel Pg^^{z) via as{Qz{l — 
z))P^^{z) with the normalization again set by probability conservation. 

Continuing in this way, we have carried out the corresponding changes for all of 
the kernels in Eq. (|6]) in the HERWIG6.5 environment, with its angle-ordered showers, 
resulting in the new MC, HERWIRI1.0(31), in which the ISR parton showers have IR- 
improvement as given by the kernels in Eq. (fT8|) . In the original release of the program. 
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V. 1.0, we stated that the time-hke parton showers had been completely IR-improved 
in a way that suggested the space-like parton showers had not yet been IR-improved 
at all. In the subsequent release, v. 1.02, the part of the space-like parton showers 
without IR-improvement associated with HERWIG6.5's space-like module HWSGQQ for 
the space-like branching process G — t- gg, a process which is not IR divergent and which 
is, in any case, a sub-dominant part of the shower, was IR-improved. In the release 
in version 1.031 the final missing IR-improvement in the space-like module hwsfbrI 
has been implemented. The IR-improvement of the module HWSGQQ in the release 
HERWIRI1.02 produces a small effect, as these considerations suggest: we see effects at 
a level comparable to the errors on the MC data in our plots when going from version 
1.0 to 1.02. In going from version 1.0 to 1.031, we do some significant effects in the soft 
Pt regime so that we recommend the use of version 1.031 in comparison with data, as we 
illustrate presently. 

We now illustrate some of the results we have obtained in comparing ISR showers in 
HERWIG6.5 and with those in HERWIRI1.0(31) at LHC and at FNAL energies, where 
some comparison with real data is also featured at the FNAL energy. Specifically, we 
compare the 2;-distributions, pT-distributions, etc., that result from the IR-improved and 
usual DGLAP-CS showers in what follows [H]. 

First, for the generic 2— )-2 hard processes at LHC energies (14 TeV) we get the com- 
parison shown Figs. |31 S] for the respective ISR 2;-distribution and distribution at 
the parton level. Here, there are no cuts placed on the MC data and we define z as 
z = -E'partonZ-E'bcam where -Ebcam IS the cms beam energy and -Eparton is the respective par- 
ton energy in the cms system. The two quantities z and for partons are of course not 
directly observable but their distributions show the softening of the IR divergence as we 
expect. 

Turning next to the similar quantities for the vr^ production in the generic 2— )-2 hard 
processes at LHC, we see in Figs. O [6] that spectra in the former are similar and spectra 
in the latter are again softer in the IR-improved case. These spectra of course would be 
subject to some "tuning" in a real experiment and we await with anticipation the outcome 
of such an effort in comparison to LHC data. 

We turn next to the luminosity process of single Z production at the LHC, where in 
Figs. [7f8|9l we show respectively the ISR parton energy fraction distribution, the Z p^ 
distribution, and the Z rapidity distribution with cuts on the acceptance as 40GeV < 
Mz, Pt > 5GeV for Z — >■ yU,"*"/!" - all lepton rapidities are included. For the energy 
fraction distribution and the pr distributions we again see softer spectra in the former 

^We thank M. Seymour and B. Webber for discussion on this point. 
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and similar spectra in the latter in the IR- improved case. For the rapidity plot, we see 
the migration of some events to the higher values of Y, which is not inconsistent with a 
softer spectrum for the IR-improved case. One might wonder why we show the Z rapidity 
here as the soft gluons which we study only have an indirect affect on it via momentum 
conservation? But, this means that the rapidity predicted by the IR-improved showers 
should be close to that predicted by the un-improved showers and we show this cross-check 
is indeed fulfilled in our plots. To understand why one has the migration to higher values 
of rapidity in the IR-improved spectra, recall the IR-improved spectra move the radiated 
partons to softer values of z and this means the produced Z's have harder values of energy 
for given pt as the px spectra are similar, and this in turn means these Z's have higher 
values of the rapidity variable. We look forward to the confrontation with experiment, 
where again we stress that in a real experiment, a certain amount of "tuning" will affect 
these results. We note for example that the difference between the spectra in Fig. [HI while 
it is interesting, is well within the range that could be tuned away by varying the amount 
of intrinsic transverse momentum of partons in the proton. The question will always be 
which set of distributions gives a better per degree of freedom. 

Finally, we turn the issue of the IR-cut-off in HERWIG6.5. In HERWIG6.5, there 
is are IR-cut-off parameters used to separate real and virtual effects and necessitated 
by the -|— function representation of the usual DGLAP-CS kernels. In HERWIRI, these 
parameters in principle can be taken arbitrarily close to zero, as the IR-improved kernels 
are integrable P^[T3]. We note that in the current version of HERWIRI, the formula for 
as{Q) is unchanged from that in HERWIG6.5 so that there is still a Landau pole therein 
and this would prevent our taking the attendant IR cut-off parameters arbitrarily close to 
zero; however, we also note that this Landau pole is spurious and a more realistic behavior 
for as{Q) as Q — )■ from either the lattice approach [IS] or from other approaches such 
as those in Refs. [ini[lZ] could be introduced in the regime where the usual formula for 
as{Q) fails and this would allow us to approach zero with the IR cut-off parameters. With 
this understanding, we now illustrate the difference in IR-cut-off response by comparing 
it for IIERWIG6.5 and HERWIRI: we change the default values of the parameters in 
HERWIG6.5 by factors of 0.7 and 1.44 as shown in the Fig. [TOl We see that the harder 
cut-off reduces the phase space significantly only for the IR-improved kernels and that 
the softer cut-off has also a small effect on the usual kernels' spectra whereas as expected 
the IR-improved kernels' spectra move significantly toward softer values as a convergent 
integral would lead one to expect; for, one must note here that the spectra all stop at 
approximately the same value zq = .00014 — .0016 which is above some of the modulated 
IR-cut-off parameters and that the peaks in the spectra are not at the respective IR 
cut-off values which are defaulted [S] at 2 = 0.000114, 0.000121 for quarks and gluons, 
respectively, as the HERWIG environment has other built-in cut-offs that prevent such 
things as the as argument's becoming too small. What the curves in Fig. [TO] show then are 
the relative "relative" probabilities for normalized spectra above zq. Specifically, the areas 
under the six curves in the figure are all equal. The curves then show the difference in 
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the shapes of the parton energy spectra for the given values of the IR cut-off parameters 
when the interplay with the HERWIG environment's other built-in cut-offs to prevent 
unphysical phenomena from occurring in the simulation is taken into account. We can 
make an estimate of the attendant relative relative probabilities as follows. We compare 
the probability P{zq < 2 < in a normalized spectrum for the spectra with the IR cut- 
offs in the usual case, where for soft gluon quanta we have dP ~ dz' /z' in usual DGLAP- 
CS theory for the IR singular part, with the analogous probability in the IR- improved 
case, where for soft gluon quanta we have dP ~ dz' j z'^~^ with 7 = 7^1, A = q, G. 
Considering the q case for definiteness, we have, for zq = 0.00015, an IR cut-off [8] 
Zcg = 0.80GeV/7TeV = 0.000114 as noted above, so that we need to compare the two 
relative probabilities 



where we also introduced notation for the HERWIG environment gluon IR cut-off |8] 
ZcG = 0.85GeV/7TeV = 0.000121 with its default value. This shows that the two rela- 
tive probabilities are in the ratio r = (1 — Zcg'^) / {—•ylia^ZcG)) so that the probability is 
enhanced in the IR-improved spectra. Putting in the numbers for we get r = 0.16 for 
the suppression of the unimproved spectrum relative to the IR-improved one. This is of 
course an overestimate of the effect because we only analyzed the IR singular terms in the 
soft regime but we see that the suppression effect is consistent with the plots in Fig. [TOl 
If we want to be more complete, we also need to analyze the q and G singular cases and 
to take the suppression of the soft gluon spectrum in the IR-improved spectrum into ac- 
count. Concerning the q case, it is the same as the q case - we get an enhancement of the 
soft region in the IR-improved spectrum relative to the unimproved one. For the gluon 
case, when the branching is G — ?■ G{z) + G{1 — z), for the part of the gluon branching 
kernel that is singular in z there is an enhancement in the relative probability that the 
gluon with fraction 1 — z will be in the soft region; the converse also holds. Finally, all 
of the soft gluon singularities are suppressed in the IR-improved cases relative to their 
unimproved cases, so that this will move particles away from the soft regime in the IR- 
improved cases relative to the unimproved cases. What is crucial to the Fig. [10] results is 
the interplay with the HERWIG6.5 environmental cut-offs such as that one on the argu- 
ment of Os = as{Qz{l — z)) where Q is the HERWIG6.5 angle-ordered evolution variable. 
This cut-off on the variable Q = z{l — z)Q means that the regime where either 2; — >■ or 
1 — 2; —J- is suppressed in both sets of spectra. But, as the unimproved spectra have their 
largest enhancements in this regime relative to the IR-improved spectra, the HERWIG6.5 
environment kills a large part of this enhancement and allows the other enhancements 
such as that in f l5^ to prevail, as we see in the figure. The behavior illustrated in Fig. [TU] 
is expected to lead to a better description of the soft radiation data at LHC and this 
expectation should be checked with experiment in the not-too-distant future. 

We finish this initial comparison discussion by turning to the data from FNAL on 




unimproved 
IR-improved 



(39) 
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the Z rapidity and spectra as reported in Refs. [IHllin]- We show these results, for 
1.96TeV cms energy, in Fig.[IIl We see that HERWIRI1.0(31) and HERWIG6.5 both give 
a reasonable overall representation of the CDF rapidity data but that HERWIRI1.031 is 
somewhat closer to the data for small values of Y . The two x^/d-o.f are 1.77 and 1.54 
for HERWIG6.5 and HERWIRI1.0(31) respectively. The data errors in Fig. [II](a) do not 
include luminosity and PDF errors |1H], so that they can only be used conditionally at 
this point. We note as well that including the NLO contributions to the hard process via 
MC@NLO/HERWIG6.510 and MC@NLO/HERWIRI1.031 [330 improves the agreement 
for both HERWIG6.5 and for HERWIRI1.031, where the xVd.o.f are changed to 1.40 
and 1.42 respectively. That they are both consistent with one another and within 10% of 
the data in the low Y region is fully consistent with what we expect given our comments 
about the errors and the generic accuracy of an NLO correction in QCD. A more precise 
discussion at the NNLO level with DGLAP-CS IR-improvement and a more complete 
discussion of the errors will appear |13]. These rapidity comparisons are then important 
cross-checks on our work. 



We also see that HERWIRI1.031 gives a better fit to the DO data compared to HER- 
WIG6.5 for low "Pt, (for < 12.5GeV, the x^/d.o.f. are ~ 2.5 and 3.3 respectively if 
we add the statistical and systematic errors), showing that the IR-improvement makes a 
better representation of QCD in the soft regime for a given fixed order in perturbation 
theory. We have also added the results of MC@NLO [33] for the two programs and we see 
that the 0[a^ correction improves the x^/d.o.f for the HERWIRI1.031 in both the soft 
and hard regimes and it improves the HERWIG6.510 x^/d.o.f for "Pt near 3.75 GeV where 
the distribution peaks. For pr < 7.5GeV the xVd.o.f for the MC@NLO/HERWIRI1.031 
is 1.5 whereas that for MC@NLO/HERWIG6.510 is worse. These results are of course 
still subject to tuning as we indicated above. 



5 Conclusions 

We have introduced a new approach to QCD parton shower MC's based on the new IR- 
improved DGLAP-CS kernels in Refs. [TUlin] and we have realized the new approach with 
the MC HERWIRI1.0(31) in the HERWIG6.5 environment. HERWIRI1.0(31) then sets 
the stage for the further implementation of the attendant [5] new approach to precision 
QED(8)QCD predictions for LHC physics by the introduction of the respective resummed 
residuals needed to improve systematically the precision tag to the 1% regime for such 
processes as single heavy gauge boson production, for example. Here, we already note 



that our new IR-improved MC, HERWIRI1.0(31), available at http://thep03.baylor.edu 



We thank S. Frixione for helpful discussions with this implementation. 
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is expected to allow for a better per degree of freedom in data analysis of high en- 
ergy hadron-hadron scattering for soft radiative effects. By comparison with the DO 
FNAL data of single Z production, we have given evidence that this is indeed the case. 
As one would expect, the integration of HERWIRI into MC@NLO is seamless, as one 
may replace HERWIG with HERWIRI directly. In both cases, MC@NLO/HERWIG and 
MC@NLO/HERWIRI, we see an improvement in the comparison with both the CDF 
rapidity data and the DO data on single Z production. We await further tests of the 
new approach, both at FNAL and at LHC 
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Figure 2: Graph of Aa{Q^,t) for the DGLAP-CS and IR-Improved DGLAP-CS kernels 
Eqs. (!26| |34|) . is a typical virtuality close to the squared scale of the hard sub-process 
- here we use = 25GeV^ for illustration. 
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Energy fraction of parton shower constituents for QCD 2^2 hard parton scattering. 
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Figure 3: The 2;-distribution(ISR parton energy fraction) shower comparison in HER- 
WIG6.5. 
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Histogram of P_T^ for QCD parton shower in I-Ierwig6.5 for 2^ 2 hard parton scattering. 
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Figure 4: The p|n-distribution (ISR parton) shower comparison in HERWIG6.5. 
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Energy fraction of n* in I-IERWIG6.5 for QCD 2^2 hard parton scattering. 
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Figure 5: The tt"^ energy fraction distribution shower comparison in HERWIG6.5. 
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Histogram of for tt^ for QCD 2^ 2 hard parton scattering. 
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Figure 6: The 7r~^ p^-distribution shower comparison in HERWIG6.5. 
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Energy fraction distribution of parton shower for single Z production. 
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Figure 7: The 2;-distribution(ISR parton energy fraction) shower comparison in HER- 
WIG6.5. 
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Generated Z transverse momentum. 
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Figure 8: The Z pT-distribution(ISR parton shower effect) comparison in HERWIG6.5. 
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The unit normalized differential cross section for Z production as a function of vector boson rapidity. 
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Figure 9: The Z rapidity-distribution(ISR parton shower) comparison in HERWIG6.5. 
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(a) 



(b) 



Energy fraction distribution of parton shower for singie Z production using DGLAP-CS kerneis. 



Energy fraction distribution of parton shower for singie Z production using IR-lmp-DGLAP-CS Iterneis 
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Figure 10: IR-cut-off sensitivity in z-distributions of the ISR parton energy fraction: (a), 
DGLAP-CS (b), IR-I DGLAP-CS - for the single Z hard sub-process in HERWIG-6.5 
environment. 
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Figure 11: Comparison with FNAL data: (a), CDF rapidity data on [Z /'-)*) pro- 
duction to e+e~ pairs, the circular dots are the data, the green(blue) hues are 
HERWIG6.510(HERWIRI1.031); (b), DO pr spectrum data on (Z/7*) production to 
e+e~ pairs, the circular dots are the data, the blue triangles are HERWIRI1.031, 
the green triangles are HERWIG6.510. In both (a) and (b) the blue squares are 
MC@NLO/HERWIRI1.031, and the green squares are MC@NLO/HERWIG6.510, where 
we use the notation (see the text) MC@NLO/X to denote the realization by MC@NLO of 
the exact 0{as) correction for event generator X. Note that these are untuned theoretical 
results. 
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